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By using the glassy (hclmm-)fluid model and boundary perturbation method, we can 
^ ■ obtain the velocity fields (as well as the flow rate; up to the second order) inside the wavy- 

I rough cylindrical nanopores which are of the same size as those samples prepared by Kim 

and Chan as well as Day et al. Our results show that the velocities measured in porous 
Vycor samples could be reproduced by carefully selecting relevant parameters but those in 
i glass capillaries are difficult to obtain. 
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1 Introduction 

(N 

> , 

0^ , A supersolid (SS) is a solid that possesses an element of superfluidity^. The existence of SS 
OO ■ ^He was first suggested 67 years ago^, while possible underlying mechanisms were discussed 30 

cn ■ 

years later"^ . The subsequent crucial step is the detection of solid He acquiring nonclassical 
. rotational inertia (NCRI)^. NCRI is indicated by a drop in the resonant period of the oscillator 
' below the transition temperature. Recently the experimental evidence for SS was reported by 
Kim and Chan in a series of torsional oscillator experiments on bulk solid helium^ where it was 
demonstrated that the decrement in the period is proportional to the SS fraction in the limit of 
low oscillation speed. 

Before the discovery of SS, there already were many researches^" about its essence and ex- 
istence. However, even after Kim and Chan announced the experimental evidence of SS, Day 
et al. conducted similar experiments^^ and showed no indication of a mass redistribution in 
the Vycor that could mimic SS decoupling and put an upper limit of about 0.003 fim/s on any 
pressure-induced SS flow in the pores of Vycor. Later on Day and Beamish's measurements^^ 
showed no indication of low temperature flow, placing stringent restrictions on SS flow in re- 
sponse to a pressure difference. The average flow speed at low temperatures was less than 1.2 x 
10~^^ m/s, corresponding to a SS velocity at least 7 orders of magnitude smaller than the critical 
velocities inferred from the torsional oscillator measurements^^"^^. Thus it was claimed that 
the origin or exact mechanism of free- flowing SS although confirmed but not yet resolved^^"^^. 
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Meanwhile, as evidenced by recent experiment (by Grigorev et al.)'^'^ where the temperature 
dependence of pressure in sohd '^He grown by the capihary blocking technique was measured; 
at temperatures below 0.3K (where the supcrsolidity was observed) they found the glassy oc 
contribution to pressure. It was then proposed by Andrccv^'^ and Rittner and Reppy^^ the 
observed supcrsolidity might be in a glassy solid (helium) state. Both idea could also be inspired 
from the annealing effect of supcrsolidity reported before^^. 

In this short paper, we offer an explanation for the argues between Day et al}'^^^'^ and Kim 
and Chan^^~^^ by treating the nanopores which were produced in porous Vycor^^ or porous 
gold samples^^ to be a kind of (cylindrical) wavy-rough nanotubes. The transport of glassy 
SS helium^^ or quantum glass^^ inside nanopores which have presumed wavy-rough wall will 
be investigated here. We adopt the verified model initiated by Cagle and Eyring^^ which was 
used to study the annealing of glass. To obtain the law of annealing of glass for explaining the 
too rapid annealing at the earliest time, because the relaxation at the beginning was steeper 
than could be explained by the bimolecular law, Cagle and Eyring tried a hyperbolic sine law 
between the shear (strain) rate : and (large) shear stress : r and obtained the close agreement 
with experimental data. This model has sound physical foundation from the thermal activa- 
tion process^^ (Eyring^^ already considered a kind of (quantum) tunneling which relates to the 
matter rearranging by surmounting a potential energy barrier) and thus it might also resolve 
the concern raised by Anderson^^ for the thermal noises to the superfiow of vortex liquid. With 
this model we can associate the (glassy) fluid with the momentum transfer between neighboring 
atomic clusters on the microscopic scale and reveals the atomic interaction in the relaxation of 
flow with dissipation (the momentum transfer depends on the activation shear volume, which is 
associated with the center distance between atoms and is proportional to RT/tq {R is the gas 
constant, T temperature in Kelvin, and tq a constant with the dimension of stress). 
To consider the more realistic but complicated boundary conditions in the wall of nanotubes, 
however, we will use the boundary perturbation technique^^ to handle the presumed wavy- 
roughness along the wall of nanotubes. The relevant boundary conditions along the wavy-rough 
surface will be prescribed below. The contents are organized into three parts : Introduction is 
firstly presented and then we describe our physical assumptions and the mathematical model. 
After that we illustrate our results and possible comparison with previous rather scattered mea- 
surements. In fact, the qualitative agreement is rather good! 

2 Physical and Mathematical Formulations 

We shall consider a steady transport of the glassy SS helium in a wavy-rough nanotube of a 
(mean-averaged radius) and the outer wall being a fixed wavy-rough surface, r = a + €sm{k9), 
where e is the amplitude of the (wavy) roughness, and the wave number : k = 2n/L. Firstly, 
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this fluid^^"^^ can be expressed as 7 = 70 sinh(r/To), where 7 is the shear rate, r is the shear 
stress, and 70 is a function of temperature with the dimension of the shear rate. In fact, the 
force balance gives the shear stress at a radius r as r = —{r dp/ dz) /2. dp/dz is the pressure 
gradient along the flow (or tube-axis : z-axis) direction. 

Introducing x = —{a/2TQ)dp/dz then we have 7 = 70 sinh(xr/a). As 7 = —du/dr [u is the 
velocity of the fluid flow in the longitudinal (2:-)direction of the nanotube), after integration, we 
obtain 

n = -Us + ^^[coshx — cosh(— )], (1) 

X a 

here, Us is the velocity over the surface of the nanotube, which is determined by the boundary 
condition. We noticed that Thompson and Troian^*^ proposed a general boundary condition for 
trasnport over a solid surface as 

Au = L07(l-^r^/^ (2) 

7c 

where Au is the velocity jump over the solid surface, is a constant slip length, jc is the critical 
shear rate at which the slip length diverges. The value of % is a function of the corrugation of 
interfacial energy. 

With the boundary condition from Thompson and Troian^^ , we can derive the velocity field and 
volume flow rate along the wavy-rough nanotube below using the verified boundary perturbation 
technique^^. The wavy boundary is prescribed as r = a -|- esm{kO) and the presumed steady 
transport is along the z-direction (nanotube-axis direction). 

Along the boundary, we have 7 = (du) / {dn)\Qj^ surface- Here, n means the normal. Let u be 
expanded in e : u = uo + eu\ + e^U2 + • ■ ■, and on the boundary, we expand u{ro + edr, 6{= Oq)) 
into 

dv"^ 

\{ro+edr,eo) = u{ro, 9) + e[dr Ur{ro, 9)] + [— 



u{r,9)\(^ro+edr,eo) = uiro,9) + e[dr Uriro,9)] + e^[^Urriro,9)] ^ = 

7«r u i./Xr 



{usUp + — [coshx - cosh(— surface (3) 



where 



^^^iiplon surface = ^Is7[(l - ^) ^^^] Ion surface' w.iipo = ^s7o[sinhx(l - (4) 

Now, on the outer wall^^ 

du ^ V(r - a — esm( k9)) . 2/,/,M-ir i 

= ^"|V(r-a-6sin(fc^))| =[1 + ^7^^°^ (^^)] ' K I (a+^<^'-,^) " 

k 

COs{k9)u0\(^a+edr,e)] = Uo^\a + e[ui^\a + Uo^^\aSm{k9)- 

k 1 k^ 

-^U0e\aCOs{k9)] + e^[--—COS^{k9)uOr.\a + U2r.\a + Ul^^\aSm{k9) + 
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-uo,,,\aSin'^ike) - ^ cos{ke){uiX + uog^\aSm{k9))] + 0{e^). (5) 

Considering L° ~ a e case, we also presume sinhx <C jc/io- With equations (1) and (5), using 
the definition of 7, we can derive the velocity field (u) up to the second order. The key point 
is to firstly obtain the slip velocity along the boundary or surface. After lengthy mathematical 
manipulations, we obtain the velocity fields (up to the second order) and then we can integrate 
them with respect to the cross-section to get the volume flow rate {Q, also up to the second 
order here) : 

rOp ra+t smikO) 

Q = / u{r, 0)rdrd9 = Q smooth + e Qpo + Qp2 ■ (6) 

JQ Jo 

Here, 9p could be 27r or 2Tr/k depending on the specific requirement for dimensionless consider- 
ation. 

3 Results and Discussions 

We firstly compare the velocity fields between the smooth and wavy-rough nanotubes in Fig. 
1 for the mean radius a = 3.5 nm (cf. Refs. 15 and 20). The wave number k is fixed to be 
10. We select the amplitude of wavy-roughness e to be 0.02 and 0.04 a (to check the geometry 
effect which is valid for small e here due to the perturbation approach). We try 70 = 10000 s~^, 
7o/7c=0.1 (cf. Ref. 29). We then set = 0.5a. The x-axis (x/a) shown in Figure 1 is the ratio 
of forcing (along the z-axis direction) per unit volume and (shear) stress. The y-axis is for the 
total velocity (u), up to the second order. Note that the real range of the (referenced) shear rate 
: 70 will depend on the specific material chosen as well as the experimental procedure^^. We can 
observe the wavy-roughness (e) will significantly enhance the flow rate^'^ due to larger surface-to- 
volume ratio and slip effect (along the wall)^^. There is net flow rate even x/a is zero (without 
forcing!)^° (it could be a persistent current or a spontaneous flow without a pressure difference 
and the velocity of this supeflow is proportional to the small amplitude of wavy-roughness, the 
(referenced) shear rate and the slip length). Our results agree with those of Ref. 21 (measured 
higher SS fraction) for larger disordered or glassy helium cases (having much larger roughness 
than that of Refs. 15 and 16). 

Once X ~ 0.1, the velocity for the smooth nanopore will be around 10 fim/s which was reported 
in Ref. 15. The upper limit of u set in Ref. 13 by Day et al. could be due to x less than 0.0001 
(presumed their sample of smooth nanoporcs, too)! Wc remind the readers that the geometry 
effect for the enhanced flow rate (for small wavy-roughness and smaller forcing) is clear as shown 
in this figure (once e is increasing) but we should keep in mind that the perturbation approach 
might limit e to be less than 0.1 a. 

To check the larger pore effect (a = 245 nm in Ref. 16), we calculate the total velocity for a=245 
nm case while keeping all other parameters almost the same. The trend (differences between 
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the smooth and rough nanotubcs) shown in Fig. 2 is almost the same as that in Figure 1. We 
can observe that, under the same selected parameters, for larger smooth nanopores, the velocity 
should be much larger than that in smaller smooth nanopores (at least 2 order of magnitudes 
larger). However, as reported in Ref. 16 for porous gold cases, the velocity is almost the same 
order of magnitude as that in Ref. 15. There might be agreement between our results and 
those in Ref. 16 for a = 245 nm ease if 70 = 10000 s^^ is lowered to less than 100 s~^! As the 
nanopore size in Ref. 14 (glass capillaries) is unknown, thus to obtain the rather small velocity 
(< 1.2 X 10~^/um/s) reported therein x should be less than 10~®. Note that those lowest curves 
illustared in Figs. 1 and 2 were obtained using the Navier-slip boundary conditions (e.g., of. 
Ref. 32). 

In brief summary, we have theoretically obtained the velocity (up to the second order) inside the 
wavy-rough cylindrical nanopores by using the glassy (solid) helium-fluid model and boundary 
perturbation method. Our results show that the calculated velocity for smaller (presumed) 
smooth nanopores (radius ~ 3.5 nm) could be of the same order of magnitude as those in Ref. 
15 after carefully selected parameters. Those measurements reported in Ref. 14, however, as 
the size of nanopores are unknown, are difficult to be reporduced using our approach! We shall 
investigate more complicated issues^^ in the future. 
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Fig. 1 Calculated velocity fields for different xl^^ (forcing (along the z-axis direction) 
per unit volume and referenced shear stress), a = 3.5 nm. e (=0.02 and 0.04 a here) is the 
amplitude of wavy roughness. The wave number of roughness (k) is 10 here. Kim and 
Chan's data^^ could be reproduced once x is around 0.1 and the slip length (L^) is 0.5 a. 
The lowest curves are only due to boundary effects (i.e., boundary slip or contributions 
from the slip length term in equation (7)) for smooth nanopores. 

As the experimental data were quite scattered and we only put the roughly averaged value 
here for comparison. 
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Fig. 2 Calculated velocity fields for different x/o (forcing (along the z-axis direction) 
per unit volume and referenced shear stress), a = 245 nm. e (=0.02 a here) is the 
amplitude of wavy roughness. The wave number of roughness (k) is 10 here. Kim and 
Chan's data^^ could be reproduced once 70 (the shear rate) is lowered down to 100 s~^. 
The lowest curves are only due to boundary effects (i.e., boundary slip or contributions 
from the slip length term in equation (7)). 



